TOK-HEP061128 



Angular Power Spectrum in Modular Invariant Inflation Model 

Mitsuo J. HayashQ 

Department of Physics, Tokai University, 
1117 Kitakaname, Hiratsuka, Kanagawa, 259-1292, Japan 

Shiro Hirafj] and Tomoyuki Takamlil 
Department of Digital Games, Osaka Electro-Communication University, 
1130-70 Kiyotaki, Shijonawate, Osaka, 575-0063, Japan 

Yusuke Okam^l and Kenji TakagQ 

Graduate School of Science and Technology, Tokai University, 
1117 Kitakaname, Hiratsuka, Kanagawa, 259-1292, Japan 

Tomoki Watanabe0 

Advanced Research Institute for Science and Engineering, Waseda University, 
3~4~1 Okubo, Shinjuku-ku, Tokyo, 169-8555, Japan 
(Dated: February 5, 2008) 

We propose a scalar potential of inflation, motivated by modular invariant supergravity, and 
compute the angular power spectra of the adiabatic density perturbations that result from this 
model. The potential consists of three scalar fields, S, Y and T, together with two free parameters. 
By fitting the parameters to cosmological data at the fixed point T = 1, we find that the potential 
behaves like the single-field potential of S, which slowly rolls down along the minimized trajectory 
in Y. We further show that the inflation predictions corresponding to this potential provide a good 
fit to the recent three-year WMAP data, e.g. the spectral index n s — 0.951. 

The TT and TE angular power spectra obtained from our model almost completely coincide with 
the corresponding results obtained from the ACDM model. We conclude that our model is considered 
to be an adequate theory of inflation that explains the present data, although the theoretical basis 
of this model should be further explicated. 
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I. INFLATIONARY COSMOLOGY 



Following three years of integration, the WMAP data 
has significantly improvedfl]. There have also been sig- 
nificant improvements in other astronomical data sets: 
analysis of galaxy clustering in the SDSS0, Q and the 
completion of the 2dFGRS|3,[3; improvements in small- 
scale CMB measurements [4 1; much larger samples of high 
redshift supernova^; and significant improvements in 
lensing data@. The constraints on cosmological param- 
eters, such as the spectral index and its running, as well 
as the ratio of the tensor to the scalar, have also been 
improved. 

The predictions of the theoretical inflation theories 
continue to find good agreement with these improved 
data sets. The ACDM model fits not only the three- 
year WMAP temperature and polarization data, but also 
small scale CMB data, light element abundances, large 
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scale structure observations, and the supernova luminos- 
ity / distance relationship [l[ • 

As a favored scenario to explain the observational data, 
it is customary to introduce a scalar field called inflaton 
into the theoretical models 0]; there are, however, several 
problems in constructing successful theories: i) Precisely 
what is inflaton? ii) What kind of theoretical frameworks 
are most appropriate to describe the theory of particle 
physics, inflation and the recently observed accelerating 
universe? iii) How can the contents of the universe be 
explained?: Baryonic matter 4%, Dark matter 23%, Dark 
energy 73% and so on. These problems seem to require a 
far richer structure of contents than that of the standard 
theory of particles. Furthermore, phenomenologically, iv) 
Is the model consistent with the observed CMB angular 
power spectra? In particular, inflaton should satisfy the 
slow-roll condition in order that the model predicts the 
nearly scale-invariant spectral index, as well predicting 
a sufficiently large number of e-folds. (See ref.p] for a 
recent review of the theories of inflation.) 

In this paper, we are proposing a potential which 
can give predictions consistent with those of the ACDM 
model. This potential was originally derived by Ferrara 
et al. in the context of T-duality and supersymmetry 
breaking in string theory via gaugino condensation. We 
have shown in ref.[Toj that inflation and supersymmetry 
breaking can occur at the same time with an appropriate 
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choice of parameters. However, because we would like to 
explore the phenomenological implications of this model 
as an inflationary theory, we will suppose in the present 
study that the potential's parameters are not restricted 
by supergravitational backgrounds. 

This paper is organized as follows: In Sec. II, starting 
with a potential form which is modular invariant in T, we 
derive a stable inflationary trajectory by fitting the pa- 
rameters. In Sec. Ill, we compute angular power spectra 
derived from this model, which are shown to agree with 
three years of WMAP data, and almost coincide with the 
corresponding results of the ACDM model. In Sec. IV, 
we present the conclusions of this study, before offering a 
discussion of these results. Finally in the Appendix, for 
completeness, we present an outline of the approach used 
to derive the present potential. In this study we will use 



Planck units such that m p i/v87r = !• 



II. INFLATIONARY TRAJECTORY AND 
STABILITY IN VARIABLE T 

In order to construct the model, we introduce the in- 
flaton field S and the gauge-singlet complex scalar fields 
Y and T, motivated by the framework of modular invari- 
ant supergravity conjectured from dimensionally reduced 
superstring s[ia [O, 111. However, we will not consider 
supergravitational backgrounds in this paper. 

We consider a scalar potential for the three fields in 
the following form, 



V(S,T,Y) = 



\Y\ 



3(S + S*)\Y\ 4 
(T + T*- \Y\ 2 ) 2 



3b 2 



T + T* - Y 
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where r\ is Dedekind's ^-function, defined by 

r](T) = e - 2rf V 24 f[ (1 - e - 27mT ), (2) 

71=1 

rj' is its derivative with respect to T, and c and b are 
the free parameters of this potential. For completeness, 
in the Appendix we will briefly review the outline of the 
derivation of Eq. ([!]) in the context of supergravity fol- 
lowing ref.|9j. In this context, b is the one-loop renor- 
malization group coefficient of gauge groups in hidden 
sectors of ref.Q, and is not free. In the present study, 
however, in order to explain observed data using this po- 
tential form, we treat b along with c as free parameters 
since we do not consider gaugino condensation. 

The potential is modular invariant in the complex 
scalar field T, and is shown to be stationary at the self- 
dual point T = I. Such T-duality often plays impor- 
tant roles in various aspects of string theories [l3j; e.g. 
this invariance is an unbroken symmetry at any order of 
string perturbation theory. One could therefore require 
the Kaler potential and superpotential to be modular in- 
variant. For simplicity, S and Y are assumed to be real 
fields and the other matter fields are neglected. 

One of the main purposes of this paper is to prove that 
the interrelation between S and Y gives rise to inflation. 
As we will see later, upon minimizing V with respect 
to Y at the fixed point T = 1, one eventually finds a 
single-field potential V(S); the scalar field S plays the 
role of the inflaton field. Usually inflaton fields must sat- 



isfy the slow-roll condition in order for the corresponding 
inflation model to be successful. Roughly speaking, c de- 
termines the energy scale of V while b determines the 
flatness of V. If c is small enough (< I), the potential 
([1]) has no local minimum in Y at fixed T other than the 
trivial minimum at Y = 0. We therefore have to choose 
the parameters c and b carefully. 

We found that the potential V(S,Y) at T = I has 
a stable minimum at (Y m i n , S m i n ) ~ (0.00877, 1.51) with 
c = 131, b = 9.4 (See Fig.l). These are the most suitable 
values for realizing the present experimental observations 
of three-year WMAP. (If we could stick to Y as a gaug- 
ino condensated scalar field, the local supersymmetry is 
broken at once with inflation, providing a seed for ob- 
servable supersymmetry breaking. The value of S m i n is 
consistent with (S + S*) = a'm^.) We can see inflation 
arises precisely due to the evolution of the scalar fields S 
and Y as follows: 

First, for the parameter values c = 131 and b = 9.4, 
the inflationary trajectory can be well approximated by 

Y min (S) ~ 0.00925e-° 0355S , (3) 

which corresponds to the trajectory of the stable mini- 
mum for both S and Y. 

In Fig. 2, we have shown a plot of V(S) minimized with 
respect to Y; for large c, this potential is asymptotically 
approximated by 

o / -S/36\ 6 

V(S) ~ - S{2S 2 + 2bS - b 2 ), (4) 

2 \ ac J 
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FIG. 1: The potential V(S, Y) at fixed T = 1 (self-dual 
point) with c = 131, b — 9.4. We can see a valley of the 
potential, Vy(S) = 0, and a stable minimum at (Ymin, S m in) ~ 
(0.00877, 1.51). 



where a is a constant determined by T: a ~ 0.824 for 
T = 1. Thus we have arrived at a single- field poten- 
tial, starting from the modular invariant potential for 
the three fields. The stability of the fixed point T = 1 
will be shown later. 
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FIG. 2: The potential V(S) minimized with respect to Y. 
The minimum value of the potential is V(Smin) ~ —5.71 x 



Next, the slow-roll parameters are defined by 

i (d a v\ 2 d a d (j v ... 

€a = 2\—) ' ^ = -^- ( 5 ) 

The slow-roll condition requires both these values to be 
less than 1. The end of the inflationary period is de- 
marked by the slow-roll parameter e a approaching the 
value 1. Beyond the end of the inflationary period, "mat- 
ter" may be produced during the oscillations around the 
minimum of the potential (reheating) at the critical den- 
sity, i.e. fi = 1. Although any successful theory of infla- 
tion should explain the mechanism of the reheating pro- 
cess, we postpone consideration of this reheating problem 
for later work. 



For the present potential the values of e$ and r]ss 
arc numerically obtained by fixing the parameters c = 
131 and b = 9.4 as shown in Fig. 3; with these parame- 
ters the slow-roll condition is satisfied. 
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FIG. 3: The evolution of the slow-roll parameters. The solid 
curve represents eg while the dashed curve denotes \r)ss\- The 
plots demonstrate that the potential V(S) is sufficiently flat. 
Inflation ends at S ~ 4.19 in our model. 



The potential V is stable at the self-dual point T = 
1 in arbitrary points in the inflationary trajectory for 
our choice of the parameters c and b. By choosing three 
points, i.e., the horizon exit, the end of inflation and 
the stable minimum, and substituting the values of S, 
Y at these points into the original V(S,Y,T), we will 
here demonstrate that the potential V(T) has a minimum 
precisely at T = 1, and hence is stable at these typical 
stages in the inflationary trajectory. The variations of 
V(T) are obtained numerically in Figs. 4 and 5 for the 
fixed parameters c = 131 and b = 9.4. 
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FIG. 4: A 3D plot around the minimum of V(T) as a function 
of the complex variable T for S m i n and Ymin. Points along 
Im T — are obviously stable. 



4 




0.5 



FIG. 5: Plots of V(T) at Im T = at three representative 
inflationary stages. It is easy to see that T = 1 is stable. 
The solid, dashed and dotted curves represent stages at the 
horizon exit, the end of inflation and the stable minimum 
respectively. 



Now we have to verify the amount of inflation. The 
number of e-folds at which a comoving scale k crosses 
the Hubble scale aH during inflation is given by0] 



N(k) ~ 62-ln 



k 1, (10 16 GeV) 4 

In 

a H 4 Vfc 



1 1 Vk ^ 
4 ln F7' (6) 



where we assume V'ond = Proh- We focus on the scale 
fc* = 0.05 Mpc~ 4 and the inflationary energy scale is 
V ~ 10~ 10 - (10 16 GeV) 4 as shown in Fig. 2. The 
number of e-folds which corresponds to our scale must 
therefore be around 57. 

On the other hand, using the slow-roll approximation 
(SRA), N is also given by 



Because n s is not equal to 1 and a s is almost negligible, 
our model suggests a tilted power law spectrum. The 
value of consistent with the recent observations; 

the best fit of three-year WMAP data using the power 
law ACDM model is n s ~ 0.951 Q. 

Finally, estimating the spectrum V-r in SRA, 



1 V s 



(12) 



we find Vn* ~ 2.36 x 10~ 9 . This result is also in 
excellent agreement with the measurements derived 
from observations. It may also be noted that the energy 
scale V ~ 10~ 10 is also within the constrained range 
obtained by Liddle and Leach [l5|. 



Gravitational waves are an inevitable consequence of 
all inflation models. The tensor perturbation (the gravi- 
tational wave) spectrum is given by 0] 



T^grav 8 



2tt J 



-V. 



(13) 



In SRA, the spectral index of V gra 
roll parameters e and r\ as 

tit = — 2e. 



3tt 2 

r is given by the slow- 
(14) 



The relative amplitude of the gravitational waves and the 
adiabatic density perturbations is given by 



N ' 



w dS - 



(7) 



We could also have obtained the number of e-folds ~ 57, 
by fixing the parameters c = 131 and b = 9.4 and inte- 
grating from 5 en d ~ 4.19 to S* ~ 11.6, i.e. our potential 
can produce a cosmologically plausible number of e-folds. 
Here S 1 * is the value corresponding to fc*. 

We can also compute the scalar spectral index and 
its running that describe the scale dependence of the 
spectrum of the primordial density perturbation V-r. = 
(H/S) 2 (H/2tt) 2 0,13; these indices are defined by 



- 1 = 



dlnVn 
dh\k 



d\nk 



(8) 
(9) 



These are approximated in the slow-roll paradigm as 

n„(S) ~ l~6e s + 2r ]ss , (10) 
a s (S) ~ 166^-244-^3), (11) 

where ^(3) is an extra slow-roll parameter that includes 
the trivial third derivative of the potential. Substituting 
S t into these equations, we have n s * ~ 0.951 and a s * ~ 
-2.50 x 10" 4 . 



(15) 



The ratio r' of the tensor quadrupole Cf rav to the scalar 
quadrupole is defined by (Peiris et al. in 0]) 



= 0.8625r = 13. 8e, 



(16) 



for the standard cold dark matter model. The gravita- 
tional wave spectrum does not evolve and remains frozen- 
in as a massless field, even after the horizon-exit, inde- 
pendent of the scalar perturbations . In contrast to 
this fact, the evolution of the primordial curvature fluc- 
tuation 1Z is given by the product of the transfer function 
T r (k) and TZ: 



?4 m) =T r {k)K k . 



Therefore, the ratio r evolves as 



V e 



Vn 



(m) 



-8T r 2 n T 



(17) 



(18) 



up to the present time. This result will be used in the 
calculation of the angular power spectra. 
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III. THE ANGULAR POWER SPECTRUM OF 
THE MODEL 

Using our model, we can calculate the angular power 
spectrum to compare with WMAP analysis and other 
experimental data[H, 0, i, H H, H, G3- The multipoles 
aim of the CMB anisotropy are defined by 



ST 



m—l 



AT = — = E a i™Yim(e), 

i>0 m=—l 

aim = / dn»AT(n)l^(e), 



(19) 



(20) 



where Yi m (n) are spherical harmonic functions evaluated 
in the direction n. The multipoles with I > 2 represent 
the intrinsic anisotropy of the CMB. If the CMB temper- 
ature fluctuation AT is Gaussian distributed, then each 
aim is an independent Gaussian deviate with 



{dir. 



= 0, 



and 



(ai m a*, m ,} = Su'Smm'Ci, 



(21) 



(22) 



where Ci is the ensemble average power spectrum, or, 
the angular power spectrum of the CMB. In general, the 
cosmological information is encoded in the standard de- 
viations and correlations of the coefficients: 



(XY) = (af m 4 m ,*) = 5 w 6 mm >Ci 



XY 



(23) 



For an arbitrary function g(x), if we use a spherical 
expansion of the form 



7(x) 



dk £ gim (k) d Z kji (kx)Y lm (0,0), (24) 

lm 



where ji is the spherical Bessel function, and 



is the 



direction of x, then the angular power spectrum Cj T and 
the temperature-polarization cross-power spectrum Cf E 
will be given by 



OF = ^fTi(k,l)Vn(k)^, 



-in 



T @ (k,l)T E (kJ)Vn(k) 



dk 



(25) 
(26) 



where Tq and Te are transfer functions and is a bright- 
ness function. 

Now we will describe the behavior of these power spec- 
tra according to our model. 

The scalar spectral index is n s (k*) = 0.951 and 
the running index is a s (fc st ) = —0.000250 at k* = 
0.05 Mpc -1 , as has already been shown. We consider a 
tensor-to-scalar ratio r' = 13.8e = 0.00910 (e = 0.000659 
at K = 0.05 Mpc" 1 ). 

We will use the CMBFAST[H|], where we have as- 
sumed the cosmological parameters to be: Qtot = 1 f° r 



the total energy density, u) = — 1 and SI a = 0.762 for 
the dark energy, Jib = 0.0418 and f2 c dm = 0.1962 for 
the baryonic and dark matter density, h — 0.73 for the 
Hubble constant. The angular TT power spectra were 
normalized with respect to 11 data points in the WMAP 
data from I = 65 to I — 210, and the same values were 
used in the analysis of the angular TE spectrum. 

By using the likelihood method[l9j. we calculated the 
X 2 values for the TT and TE spectra, and also for their 
total sum. The results are shown in Table I. 



The x 2 values for the ACDM model with t = 0.1 
which were given by Spergel et al. [l[ , were also calculated 
by the same method. (For the one-year WMAP data, 
t = 0.17 was favored.) On the other hand, the best fit of 
our model is realized at r = 0.090 for both TT and TE 
modes, which falls within experimental error, while the 
X 2 values of our model seem to be better than those of 
the ACDM model. 

The angular power spectrum of our model for the TT 
mode at r = 0.090 are presented in Fig. 6. Figs. 7 and 8 
show the TE spectra for I < 50 and with more detailed 
data I < 20 respectively also at t — 0.090. Note that 
both the TT and TE spectra almost completely coincide 
with those of the ACDM model across the whole range 
of I. We would like to emphasize these two results as 
characteristic features of our model; within the area of 
scalar inflation models, our model can be regarded as an 
alternative to the ACDM model. 

Although both our model and the ACDM model per- 
form on the whole satisfactorily in explaining the WMAP 
data, there remains one inconsistency in the suppression 
of the spectrum at large angular scales (I = 2)P, [l7T |. 
This problem is at present left to future investigations. 

In summary, the model we have here investigated is 
consistent with the present observational data and the 
ACDM model. 



TABLE I: The x 2 values for the TT and TE spectrum and 
their total sum. The values of the ACDM model were calcu- 
lated by using the results of Spergel et aZ.[l[. 





T 


TT 


TE 


Total 


Our model 


0.090 


1057.17 


418.50 


1475.67 


ACDM 


0.088 


1057.56 


418.50 


1476.06 
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FIG. 6: Temperature angular power spectrum (TT). 
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FIG. 7: Temperature-polarization cross-power spectrum 
(TE) for I < 50. 
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FIG. 8: Temperature-polarization cross-power spectrum 
(TE) for I < 20 with more detailed data. For larger val- 
ues of I the results of our model are almost indistinguishable 
from those of the ACDM model. 



IV. CONCLUSION AND DISCUSSION 

We conclude that the results obtained with the present 
model are in good qualitative and quantitative agreement 
with both the angular power spectra of the three-year 
WMAP data, and the corresponding results obtained 
from the ACDM model. 

Although our model started from considering a poten- 
tial of three fields, it could eventually be treated as a 
potential of a single field. It is interesting to speculate 
whether the present model may be regarded as a model 
of multi-field inflation, including the effects of fluctua- 
tions in the three fields. From a particle physics point 
of view, it seems more natural to expect more than one 
scalar field to roll during inflation. In this case it may be 
necessary to consider a spectrum of isocurvature as well 
as curvature, and the correlations between the two(20|. 
Thus our model may be seen as a prelude to a full multi- 
field inflation model. 

As we mentioned before, the inflaton potential we 
used was originally derived by Ferrara et al. 0] in the 
context of dimensionally reduced supergravity from su- 
perstring theories. Their potential form was based on 
the idea of constructing an effective theory of gaugino 
condensation [2 1| . incorporating the target-space duality, 
where the gaugino condensation has been described by 
a duality-invariant effective action for the gauge-singlet 
gaugino bound states Y coupled to the fundamental 
fields as the dilaton S and moduli T. On the other 
hand, in ref.[22j], the gaugino-condensate has been re- 
placed by its vacuum expectation value to yield a duality- 
invariant "truncated" action that depends on the funda- 
mental fields only. The equivalence between these two 
approaches has been proved in refs. [23[ . 

It appears that supergravity is one of the most plau- 
sible frameworks to explain new physics, including un- 
detected objects, such as the inflaton, dark matter and 
dark energy. In particular, since the inflaton field is con- 
cerned with Planck scale physics, a dilaton field seems to 
be the most likely candidate for the inflaton. The con- 
struction of a realistic supergravity is another problem 
that must be tackled in the future, and would appear to 
be a fruitful approach. 

In the present study we have treated the parameters b 
and c as completely free parameters. However, if we hope 
to consider our model as one given by string-inspired su- 
pergravity with gaugino condensation of the hidden sec- 
tor, the value b = 9.4 is too large to be realistic, be- 
cause b is defined by gfcy as the coefficient of the one- 
loop /3-function of the renormalization group equation 
for a gauge coupling constant of a gauge group, e.g. Eg. 
Therefore our conclusion is restricted to presenting our 
model as a possible form of potential which gives rise to 
an adequate inflation consistent with the present WMAP 
data. 

Because the agreement with the WMAP observations 
does not seem merely accidental, our next tasks include 
the investigation of an alternative derivation of an infla- 
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tion potential of similar form to the present model using 
supergravitational theory. Furthermore, the reheating 
and matter production (dark and baryonic matters) fol- 
lowing inflation will also be an immediate further area of 
investigation^, 24]. 



APPENDIX: DERIVATION OF THE SCALAR 
POTENTIAL 

In this appendix we briefly review the derivation of 
the scalar potential Eq. |T]) following Ferrara et aL[9|] in 
the context of modular invariant supergravity. Assuming 
that the compactification of the superstring theory pre- 
serves N = 1 supersymmetry, an effective theory should 
be of the general type of N = 1 supergravity coupled to 
gauge and matter fields. The most general form of the 
Lagrangian in N = 1 supergravity at the tree-level is [ill]: 



e- K / 3 S S 



-[siw] F +[f ab w^w^] F , 

( A -l) 

where the Kahler potential K is given by 



K = -In (S + S*) - 3 In (T + T* 
and the gauge function / a & is 

fab — 8 a bS- 



1$,: 



(A.2) 



(A.3) 



In order to construct an effective theory of gaugino con- 
densation, we introduce the composite superfield Y of 



the gaugino condensation [9t |21|: 

Y 3 = 5 ab W^Wl/Sl = (AA + ■ • • )/S 3 , 



(A.4) 



where A is the gaugino field in the Hidden sector. 

The effective Kahler potential and superpotential in- 
corporating modular invariant one-loop corrections are 
given by[9( 

K = -ln(S + S*)-3ln(T + T* - \Y\ 2 - ^.f) , (A.5) 
and 



W = 3bY 3 In 



ce s/3b r7?2(r) 



W n 



(A.6) 



where r\ is Dedekind's 77-function, c is a free parameter in 
the theory and b = g|pj (/?o is the one-loop beta-function 



coefficient). Since (S + S*) = a'm^u the choice 



\e- K/3 S So 



8=6=0 



[S + S] t 



=8=0j 



(A.7) 



corresponds to the conventional normalization of the 
gravitational action: 



•Cgrav ~ [e K ^ 3 SoSo]g =e=0 R. 



(A. 



Then, after substituting Eqs. (|A.5|) - (|A.7[> into Eq. (|A.l|) 
and eliminating the auxiliary fields, the scalar potential 
is obtained as Eq.JT]). 
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